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Abstract We report on the experimental observation of vortex formation and production of 
tangled vortex distribution in an atomic BEC of *'^Rb atoms submitted to an external oscil- 
latory perturbation. The oscillatory perturbations start by exciting quadmpolar and scissors 
modes of the condensate. Then regular vortices are observed finally evolving to a vortex tan- 
gle configuration. The vortex tangle is a signature of the presence of a turbulent regime in the 
cloud. We also show that this turbulent cloud has suppression of the aspect ratio inversion 
typically observed in quantum degenerate bosonic gases during free expansion. 
PACS numbers: 03.75.Lm, 03.75.Kk, 47.37.-Hq, 47.27.Cn 



1 Introduction 

Superfluidity and the properties associated with it like quantized vortices and turbulence, 
or Quantum Turbulence (QT), have been extensively studied in superfluid Helium below 
the Lambda point both experimentally and theoretically since it has been discovered about 
50 years agoi. QT is a phenomena characterized by the appearance of quantized vortices 
distributed in a tangled way^'^'^. Until recently, turbulence in He-superfluid could only be 
observed by indirect methods^, since individual vortices are too small to be observed di- 
rectly. A few years ago QT in Helium* has been directly observed by means of imaging 
solid particles or atoms trapped in the vortices cores. Theoretically, QT in Helium has been 
satisfactorily modeled by the vortex filament model ^, though some features such as vor- 
tex reconnections must be introduced artificially. Up to very recently. Helium was the only 
superfluid providing evidences of QT and revealing its properties. 
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The achievement of BEC in trapped atomic gases^ and the subsequent observation of 
quantized vortices in these samples^ opened up the possibility to study turbulence in a more 
controlled fashion, since one can control the main characteristics (interaction, atomic den- 
sity, number of atoms, trapping configuration) of the atomic cloud. In fact, the study of QT 
in atomic superfluids may shine light on the turbulence characteristics that are due to the su- 
perfluidity as a macroscopic quantum state and reveal which of them are in close connection 
to the superfluidity itself. 

Quantized vortices in BECs were preferentially produced by inducing rotation in the su- 
perfluid atomic cloud around a single axis by means of rotating asymmetric trapsiS and/or 
stirring potentials ' ' . That methods have little chances to produce QT because the single 
sense of rotation and consequent regularity with which the vortices are nucleated. In fact, it 
has been shown experimentally— and theoretically— that rotating a BEC in a single direc- 
tion generate vortices that evolve to the well-known vortex lattices. Recently, vortex forma- 
tion has been also observed by merging multiple trapped BECs and by coherent transfer 
of orbital angular momentum from optical fields to the condensed sampl&i^. 

In this communication, we report on the production of multidirectional vortices which 
allowed the experimental observation of QT in a magnetically trapped BEC of '^^Rb atoms 
evidenced by the observation of tangled vortices in the quantum sample. Following the iden- 
tification of QT, a dramatic change in the hydrodynamic behavior of the atomic cloud during 
free expansion was revealed. The usual aspect ratio inversion of the condensate is suppressed 
and the turbulent cloud expands keeping its anisotropic initial spatial profile. We strongly 
believe that this behavior is intrinsically associated to the presence of turbulence in the super- 
fluid. That effect resembles a kind of self-trapping due to the tangled vortices configuration. 



2 Experimental Procedure 

The description of the experimental sequence to achieve Bose-Einstein Condensation as weU 
as the protocol to generate vortices in the condensate are described in details in Refs— ii2ili. 
In brief, we produce a BEC of ^^Rb containing (1—2) x 10^ atoms in a cigar-shaped mag- 
netic trap with frequencies given by CO,- = 2k x 210(3) Hz and a>x = 2k x 23(3) Hz. We 
measure our trapping frequencies using the standard procedure with cold atoms. We induce 
simple dipolar motion of the condensate by rapidly changing the bias field of the trapping 
potential, let the condensate evolve freely and measure its position as a function of time. 
That movement is periodic and its frequency is the trapping frequency. Recently, we have 
published a theoretical analysis of the expansion dynamics of our condensate, when it co- 
exists with a large thermal component, and with those above mentioned trapping frequencies 
as input for the model, the experimental data can be quite well modeled. FiglTJa) illustrates 
the experimental setup, showing the trapping coils and absorption image beam. 

After reaching BEC, an extra field, produced by a pair of anti-Helmholtz coils is su- 
perimposed on the magnetic trap field as shown in Fig[TJb). The axis of the extra coils is 
close to the trap axis (the angle between these axes is 5 ~ 5°). The center of the extra field, 
defined by the zero-field amplitude position, is close to the QUIC trap minimum. An oscil- 
latory current is applied to the coils, always having the same sign and always starting from 
zero, so we do not give an abrupt kick to the condensate in the beginning of the excitation 
phase (/(.„,■/ = Iq [1 — cos(i2rj]). The anti-Helmholtz coils produce a quadrupolar magnetic 
field, given by B = A (— 2jcx' +yy' +zz'), where A is the magnetic field gradient and x' is 
the direction parallel to the line that join the center of the coils. 
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Fig. 1 Schematics of (a) above view, and (b) side view of the imaging beam and the coil aiTangement in our 
system, showing the QUIC trap coils and the excitation coils (AC coil). The drawing is not to scale and the 
excitation coils misalignment is made bigger for clarity. 

The experimental sequence is then as follows: finished the evaporative cooling, the ex- 
ternal field is turned on from zero to 100 ms. For this communication the frequency of 
oscillation has been kept constant at 200 Hz and the amplitude of the magnetic field gradi- 
ent has been varied from zero up to 190 mG/cm. After the end of the oscillation stage, atoms 
are left trapped for extra 20 ms before being released and observed in free expansion by a 
standard absorption imaging technique. 



3 Results 

3.1 Vortex nucleation and emergence of turbulence 

As the oscillations are introduced in the cloud, for small amplitudes of the oscillating field 
as well as short excitation periods we observe dipolar, quadrupolar and scissors modes of 
the BEC. The importance of this experimental observation is that the excitation of these 
modes give evidence that our extra field disturbs mechanically the cloud, changing not only 
its rest position but also its symmetry axis and shape by means of changing the trapping 
frequencies. The condensed and thermal clouds do not follow these mechanical disturbances 
in the same way giving rise to a relative movement between them. That relative movement is 
a fundamental feature for the possible mechanism of vortex formation that we briefly discuss 
below. Increasing both parameters (time of excitation and amplitude) we start to observe the 
formation of vortices. They seem to be formed preferably on the edges of the quantum fluid. 
The number of vortices observed grows as a function of the amplitude and/or duration of 
excitation. Examples of such regimes are shown in Fig|2] The formation dynamics of the 
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vortices as a function of the amplitude of the excitation field has been extensively discussed 
in Ref.'^. 

When the amplitude and/or excitation time is increased above 160 mG/cm a completely 
different regime takes place. We observe an increase in the number of vortices followed by 
a proliferation of vortex lines not only in the original plane of the individual vortices, but 
distributed for the whole sample, covering many directions. An example of such a regime 
can be seen in the images of the right column of Fig[3^. 

The presence of quantized vortices non regularly distributed along the whole sample in 
all directions is quite characteristic of turbulence in the quantum fluid and it is here taken 
as evidence for this regime. In fact, our way of nucleating vortices, where non-equivalent 
translational and torsional movements of the quantum cloud occur in different planes of 
symmetry of the sample, can be seen as analogous to what is proposed in/* for combined 
rotations. 

Compared to liquid Helium, the total number of vortices contained in the cloud is quite 
limited. The relative low density of the atomic fluid together with its finite size cloud creates 
a limit for the maximum obtained number of vortices that is much smaller than in liquid 
Helium. Consequently some turbulence features might not be possible to be observed and/or 
be very different compared from that obtained in liquid Helium. In this class are effects 
related to large volumes cases of quantum turbulence, containing many vortices, specially 
those related to the Komolgorov spectra, where large volumes are needed. On the other hand, 
the low vortex density may allow one to make detailed observations of vortex reconnections 
and all its related phenomena such as cascade-like processes. In order to test the Komolgorov 
spectra, the Bragg scattering technique used by Muniz et al— to study vortex lattices is 
probably the best candidate. In the case of observation of turbulence decay, in-situ, non- 
destructive images should be used^'. 

Concerning the vortex formation during oscillating perturbations, we believe that a pos- 
sible mechanism of vortex formation are the so-called Kelvin-Helmholtz instabilitieSiS, a 
well-known phenomena in the scope of fluid interface theory and experiments. These insta- 
bilities occur at the interface of two fluids with a relative motion and give rise to the for- 
mation of vortices. This phenomenon has been theoretically investigated^-' in normal fluid- 
superfluid mixtures and directly observed in mixtures of A and B phases of superfluid liquid 
Helium^"*. Nevertheless, it has never been observed in atomic BEC, though we have shown 
strong evidence of it in a previous work'^. In our case, the superfluid is the BEC compo- 
nent whereas the normal fluid is the thermal cloud. Both clouds have different frequencies 
of oscillation for quadrupolar and scissors modes, giving rise to a relative movement and 
consequent friction between them. Since the relative movement is periodic, we believe that 
at this point vortices and anti-vortices are nucleated, leading the cloud through the regular 
vortices (with possible presence of anti-vortices) regime until it reaches the vortex tangle 
regime when this number increases dramatically. The presence of anti-vortices is evidenced 
by the evolution to turbulence, which is not possible without a mixture of both types of 
vortices. In fact, if only one type of vortex were present, the cloud would evolve to a vor- 
tex lattice. Also, we observe formation of vortex clusters^ similar to what is theoretically 
predicted^^ for the coexistence of vortices and anti-vortices in the sample. 

3.2 Unusual free expansion 

The appearance of tangled vortices is certainly a strong evidence of turbulence in our system. 
Nevertheless, as expected from any system that makes the transition from a regular flow to 
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Fig. 2 (Color online) Top line: evidences of scissors modes (bending of the symmetry axis) after low- 
amplitude excitation by the external magnetic field. Bottom line: from left to right, increasing number of 
vortices observed in the BEC as the amplitude of the excitations is increased. All images taken after 15 ms of 
free expansion. 



a turbulent flow, the hydrodynamic behavior of the cloud is also changed. That unusual 
behavior manifests itself as a suppression of the aspect ratio (ratio between the most and 
least confining axis of the cloud) inversion of the cloud. This can be seen clearly in Figs[3t 
and[3l3 where we place side by side an expanding non-excited BEC where the aspect ratio 
inversion is seen and a turbulent cloud that expands keeping its aspect ratio constant. FiglSj) 
is the numerical evaluation of the aspect ratio of the expanding clouds as a function of time. 
It is important to emphasize that a non-condensed cloud expands towards isotropy or, in 
other words, it starts with an anisotropic shape (the shape of the confining potential), but it 
expands and reaches unitary aspect ratio, does not going beyond. Those facts together let us 
conclude that the turbulent cloud keeps a quantum behavior while expanding, though this 
behavior is not the expected one for a quantum degenerate bosonic cloud. 

In brief, the equation that governs the behavior of a BEC is the Gross-Pitaevskii equa- 
tion, given by 



2m 



+ y,,,(r) + «7o|f(r,Ol' 



f(r,f) = iH^f(r,f), (1) 



where —^5^ is the kinetic energy, Vf..v/(r) is the confining potential and {/o|f (r,f)|^ ac- 
counts for the interaction energy. In a expanding condensate, V'e.tf(r) = and the whole 
expansion is governed by the interplay of kinetic and interaction energies. For a thermal 
cloud the kinetic term dominates and since it is isotropic the expansion is also isotropic. 
For a BEC we have the other way round behavior: —^5^ << t/o |'f'(r,f)l^' which means 
that the anisotropic interaction energy (coming from the anisotropic spatial distribution 
|f (r,f)|^) governs the expansion, giving rise to the aspect ratio inversion phenomena, one 
of the strongest signatures of quantum degeneracy in a system of bosons. In Ref.^^ authors 
present a detailed description of this problem for a J = bosonic cloud. For T > we re- 
fer to Ref. where a modified Hartree-Fock treatment takes into account the effects of the 
thermal cloud on the expanding condensate. In any case, the aspect ration inversion of the 
condensate is predicted. Under these considerations, we cannot explain the observed behav- 
ior of the turbulent cloud. There is certainly a kind of self-trapping due to the vortex tangle 
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Fig. 3 (Color online) (a) Side-by-side images of a regular, non-excited BEC and a turbulent cloud showing 
the aspect ratio inversion on the former case and the suppression of that inversion in the turbulent regime, (b) 
Aspect ratio (ration between main axis) of the BEC and turbulent clouds evidencing the inversion of the first 
and maintenance of the latter. 

configuration that prevents the cloud to leave its initial spatial distribution. It may well be 
that the presence of a randomic distributed field of rotations prevents the atoms to turning 
into an isotropic expansion. The idea is that the rotational field distribution creates an ef- 
fective potential that expands together with the cloud. We should also observe that the time 
of flight for the turbulent cloud indicates that a much hotter fluid is obtained. That can be 
an indicative from decay mechanisms of the QT regime. Formation of Kelvin waves during 
vortices reconnections may be present and producing energy release in the form of phonons 
that increase the overall cloud energy. The increase in temperature of the cloud is a topic 
currently under investigation. 



4 Conclusions 

In this work we have described a technique no generate vortices and anti-vortices in a BEC 
with evolution to Quantum Turbulence. Turbulence manifests itself by the appearance of a 
vortex tangle in the sample in contrast to the regular vortices regime. The turbulent state 
shows up also in a dramatic change in the hydrodynamic behavior of the cloud, where a 
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suppression of the aspect ratio inversion is observed. Instead, the cloud expands keeping its 
original spatial profile and constant aspect ratio. 
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